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Abstract 

We study two-dimensional cellular automata, each cell takes three states: rest- 
ing, excited and refractory. A resting cell excites if number of excited neighbours 
lies in a certain interval (excitation interval). An excited cell become refractory 
independently on states of its neighbours. A refractory cell returns to a resting 
state only if the number of excited neighbours belong to recovery interval. The 
model is an excitable cellular automaton abstraction of a spatially extended semi- 
memristive medium where a cell's resting state symbolises low-resistance and 
refractory state high-resistance. The medium is 5emz-memristive because only 
transition from high- to low-resistance is controlled by density of local excita- 
tion. We present phenomenological classification of the automata behaviour for 
all possible excitation intervals and recovery intervals. We describe eleven classes 
of cellular automata with retained refractoriness based on criteria of space-filling 
ratio, morphological and generative diversity, and types of travelling localisations. 

Keywords: cellular automaton, excitable medium, space-time dynamics 



1. Introduction 

The memristor (a passive resistor with memory) is a device whose resistance 
changes depending on the polarity and magnitude of a voltage applied to the de- 
vice's terminals and the duration of this voltage's application. Its existence was 
theoretically postulated by Leon Chua in 1971 based on symmetry in integral 
variations of Ohms laws ||71|8l|9l|. The memristor is characterised by a non-linear 
relationship between the charge and the flux; this relationship can be generalised 
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to any two-terminal device in which resistance depends on the internal state of 
the system [[8J. The memristor cannot be implemented using the three other pas- 
sive circuit elements — resistor, capacitor and inductor - therefore the memristor 
is an atomic element of electronic circuitry [|71[8l|9l. Using memristors one can 
achieve circuit functionalities that it is not possible to establish with resistors, ca- 
pacitors and inductors, therefore the memristor is of great pragmatic usefulness. 
The first experimental prototypes of memristors are reported in [fT8l[T0ll20l . Po- 
tential unique applications of memristors are in spintronic devices, ultra-dense 
information storage, neuromorphic circuits, and programmable electronics fn\. 

Despite explosive growth of results in memristor studies there is still a few 
(if any!) findings on phenomenology of spatially extended non-linear media with 
hundreds of thousands of locally connected memristors. We attempt to fill the gap 
and develop a minimalistic model of a discrete memristive medium. The only — 
so far — approaches to develop cellular automata model of a memristive medium 
are Itoh-Chua memristor cellular automata, where cellular automaton lattice is 
actually designed of memristors [fT3l , and Adamatzky-Chua model of memrisi- 
tive cellular automata based on structurally-dynamic cellular automata |l6l. Both 
models imitate memristive properties of links, connections between cells of au- 
tomata arrays but not the cells themselves. In this paper we explore the scenario 
when link between cells are always 'conductive' but cells themselves can take 
non-conductive (refractory) states. 

We define semi-memristive automata, or excitable cellular automata with re- 
tained refractoriness, in Sect.[2j Section 2.5 presents grouping of cell-state transi- 
tion functions into eleven classes. Hierarchies of classes based on their morpho- 
logical diversity, space-filling ratio and expressiveness are constructed in Sect. [3] 
In Sect. |4] we overview a 'zoo' of travelling localizations observed in simulated 
cellular automata and exemplify interactions between the localizations. Overview 
of the results is provided in Sect. [5j 



2. Definitions and methods 

2.1. Semi-memristive automaton 

A cellular automaton is an orthogonal array of uniform finite-state ma- 
chines, or cells. Each cell takes the finite number of states and updates its states 
in discrete time depending on states of its closest neighbours. All cells update 
their states simultaneously by the same rule. We consider eight-cell neighbour- 
hood and three cell-states: resting o, excited +, and refractory — . Let u{x) = {y : 
|x — jIloo = 1} be a neighbourhood of cell x. Let = 'Lyeu{x)X{y' i+) be the sum 
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of excited neighbours of cell x at time step t, where +) = 1 if / = + and 

X{y' ,+) = otherwise. Let & = [9i, 62], I < 0\ < 62 < 8, be an excitation inter- 
val, and 4> = [^1,^2], 1 < 01 < 02 < 8, be a recovery interval. A cell jc updates its 
state by the following rule: 



A resting cell becomes excited if the number of its excited neighbours lies in the 
interval 0. An excited cell takes refractory state unconditionally. A refractory cell 
returns to its resting state if number of its excited neighbours lie in the interval 4>. 

In computer simulations we consider 0i = 1 or 2 and 0i < 02 < 8; there is 
no point to consider 0i > 2 because for such lower boundary of excitation inter- 
val most excitation patterns quickly extinct. We adopt the same boundaries of 
recovery interval 4>: 0i = 1,2 and 0i < 02 < 8. Further we refer to excitation and 
retention of refractoriness functions as ^{Oi , 02, 01 5 02)- 

2.2. Why semi-memristive? 

Why do we call the automaton described above a semi-memristive automa- 
ton? Because we assume that the refractory state symbohses high-resistivity, or 
non-conductivity. A refractory cell can not be excited whatever number of excited 
neighbours it has, therefore it can not 'conduct' excitation. In a fully memris- 
tive automaton both transitions from excited to refractory and from refractory to 
resting must be controlled by a local excitation density. In the semi-memristive 
automaton, the transition from excited to refractory is unconditional and happens 
independently of a local excitation density but the transition from a refractory state 
to resting state depends on a local excitation density. 

2.3. Experiments 

We experiment with cellular automaton array of nxn cells, n = 500, and 
absorbing boundary conditions. While testing automata's response to external 
stimulation we use point- wise and D- stimulation (a spatially extended stimulation 
in a disc). By point- wise stimulation we mean excitation of one cell (for func- 
tions, where 0i = 1) or two neighbouring cells (for functions, where Gi — 2) of 
elsewhere resting cellular array. The D-stimulation is implemented as follows. 
Let D-disc be a set of cells which he at distance not more than r from the centre 
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(n/2,n/2) of the cellular array. When undertaking D-stimulation we assign ex- 
cited state to a cell of D with probability 0.1. In computer experiments we used 
r = 100. 

2.4. Characteristics used for classification 

To classify morphology of configurations generated by functions ^{Q\ , ^, 0i , 02), 
01 = 1,2, 01 < 02 < 8, 01 < 02 < 8, we D-stimulate cellular automaton lattices, 
and subdivided them into classes by visually detected likeness of excitation and 
refractory patterns generated and based on the following characteristics: 

• Space-filling ratio. Let automaton of n x n cells being D-stimulated with 
disc radius r, then space-filling ratio is a ratio of non-resting cells in con- 
figuration recorded at time step n — r (i.e. when front of propagating distur- 
bance just reaches edge of the automaton array) to the size of array n^. 

• Morphological diversity is the number of proportion of different neigh- 
bourhood configurations. It is evaluated using Shannon entropy and Simp- 
son's index. Let W = {o, -|-, — } be a set of all possible configurations 
of a 9-cell neighbourhood w{x) = u{x) Ux, x G L. Let c be a configura- 
tion of automaton, we calculate number of non-resting configurations as 
rj = Y.xeL^i^)y where e{x) =0 if for every resting x all its neighbours 
are resting, and e{x) — I otherwise. The Shannon entropy is calculated as 
— LweVF (^(^) /'H ■ ln{v{w) /t])), where v(w) is a number of times the neigh- 
bourhood configuration w is found in automaton configuration c. Simpson's 
index is calculated as 1 — Lweiv(^(^) I^Y- 

• Generative diversity is a diversity of a configuration developed from a point- 
wise excitation. 

• Expressiveness of a function is calculated as the Shannon entropy divided 
by space-filUng ratio. Generative expressiveness is the expressiveness cal- 
culated after point-wise stimulation. 

• Power of a class is the number of functions in this class. 

• Ratios of excited and refractory states are calculated as n^^Y^xeiXi^ i'^) 
and n^^Y.xei.X{^' ■) ")> where xi^-i^) = 1 if « = ^, and xi^-i^) — other- 
wise. 
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2.5. Classes 

We separate functions 62, ^1,^2) into eleven classes, characterise mor- 

phology of configurations, generated by automata from each class, and provide 
details on membership of each class. 

2.6. Class Ci 

D-stimulation gives birth to propagating quasi-chaotic pattern of excited and 
refractory states, the pattern fills the automaton array (Fig. [1^). A front of the 
pattern is comprised of branching localised excitations, which periodically merge 
in a single excitation wave-front; the pattern is almost isotropic with slight in- 
dentations along north-south and west-east axis. Excited and refractory states 
are exhibited equally in patterns generated by automata on class Ci (Fig. |2} Ci). 
The class Ci is comprised of 27 functions: ^{\\\a), a= 1,...,8; ^(112a), 
a = 3,...,8;^(121a),a = 2,...8;^(122a),a = 3,...8. 

2. 7. Class C2 

D-stimulation leads to propagating quasi-chaotic pattern with domination of 
refractory states. The pattern fills the automaton array (Fig. [I]?). The pattern's 
propagating front is comprised of travelling localizations, which branch periodi- 
cally. Excited localizations travelling along north-west-south-east and north-east- 
south-west axis leave a distinctive trail of refractory states (Fig. [I]?). There is a 
tail of gradually extinct excitations. These spatially extended excitations collapse 
into localised oscillating localizations, most of which extinguish with time (Fig.[2| 
C2). The class C2 has three functions ^(1122), ^(1211), ^(1222). 

2.8. Class C^ 

Initial random stimulation leads to formation of a single circular excitation 
wave, followed by trains of excitation waves travelling north-west, north-east, 
south-west and south-east. The wave-trains are followed by discoidal growing 
quasi-random pattern of excited and refractory states (Fig. [1]:). North, south, 
east and west domains of cellular array lying between the growing patterns are 
filled entirely with refractory states. These refractory domains contribute towards 
slight prevalence of refractory states (Fig. [2} C3). The class C3 has eight functions 
^ilalb), a = 3, . . . , 6, = 1,2. 
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(a) Ci 



(b) C2 




(c) C3 (d) C4 



Figure 1: Exemplars of configurations generated by functions from classes Ci ... C4 with D- 
stimulation. Excited cells are red (light grey), refractory cells are blue (dark grey) and resting cells 
are white. 
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Ratio of exdted cells 



Figure 2; Ratio of excited states to resting states vs ratio of refractory states to resting states in 
configurations developed by a cellular automat with retained refractoriness with D-stimulation. 
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2.9. Class C4 

The D- stimulated automata response is characterised by a single circular ex- 
citation wave, which encapsulates wave trains propagating north-west, north-east, 
south-west and south-east, and envelopes of wave-fronts travelling north, east, 
south and west (Fig.[T|l)). The wave envelopes in C4 occupying the same space 
domains are refractory domains in C3. This is why a cluster of C4 functions is 
positioned symmetrically to cluster of C3 functions with respect to the diagonal of 
equal ratios of excited and refractory states in Fig. [2} The class C4 has 48 functions 
^{\abc), a = 3, . . . ,6, b = 1,2, c = 3,...,8. 

2.10. Classes C5 and C(, 

In response D-stimulation automata from these classes generate circular waves 
of excitaton (Fig.[3^b). The domain D of initial perturbation is a random pattern of 
refractory and resting states. No excitation is observed inside D in automata from 
class Cg and usually just few localised oscillating excitations in automata from C5. 
The circular wave-front leaves behind an almost uniform field of refractory cells 
with radial traces of resting states towards north-west, north-east, south-west and 
south-east directions (Fig.|3^b). Functions from classes C5 and Cg are positioned 
close to class C2 in Fig. |2] due to high ratio of refractory states in configurations 
they generate. The class C5 has four functions ^(la22), a = 3, ... ,6 and the 
class C(, has six functions ^(labc), a = 1,S, b = 1,2, b < c <2. 

2.11. Classes C-j and Cg 

Automata from these classes generate a single circular excitation wave while 
perturbed by D and a domain of refractory and resting states inside boundaries 
of D. In automata from Cg a propagating excitation wave leaves a trail of refrac- 
tory states along north-west- south-east and north-east-south-west axis, the rest of 
cells remain in the resting state. Propagating wave-front leaves no traces in au- 
tomata from Cg (Fig.[3]:d). Class C7 has 16 functions ^(lalZ?), a = 7, 8, = 3,4; 
J^{\a2b), a = l,%, = 3,...,8; and, class Cg has eight functions ^{lalb), 
a = l ,% and = 5, . . . ,8. 

2.12. Class Cg 

D-stimulation of an automaton from Cg produces a domain of refractory and 
resting states (inside boundaries of D) and localizations travelling along north- 
south and west-east axis (Fig. |4^). The localizations leave traces of refractory 
states. Configurations produced by functions from Cg has the lowest (amongst all 
other classes) ratio of excited states (Fig. [2} Cg). The class Cg has 17 functions 
J^{22ab), a=\,2,a<b<%\ ^(2all), a = 3,4. 
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(a) Cs 



(b) Cfi 




(c) Cy (d) Cg 



Figure 3: Exemplars of configurations generated by functions from classes C5 . . . Cg with D- 
stimulation. Excited cells are red (light grey), refractory cells are blue (dark grey) and resting cells 
are white. 
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(a) C9 



(b) Cio 




(c) Cn 

Figure 4: Exemplars of configurations generated by functions from classes Cg ... Cn with D- 
stimulation. Excited cells are red (light grey), refractory cells are blue (dark grey) and resting cells 
are white. 
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Class 



Shannon entropy Simpson's index Space-filling ratio Expressiveness 



Ci 

C2 

C3 
C4 
C5 

C6 
C7 

Cg 

C9 

Clo 
Cii 



7.01 
5.06 
5.99 
5.55 
1.79 
2.15 
3.95 
4.56 
3.89 
5.20 
4.27 



1.0 
0.95 
0.96 
0.98 
0.53 
0.65 
0.96 
0.99 
0.95 
0.98 
0.97 



0.94 
0.95 
0.96 
0.96 
0.96 
0.96 
0.42 
0.15 
0.14 
0.60 
0.27 



7.41 
5.35 
6.24 

5.8 

1.87 
2.23 

9.5 
31.21 
27.69 

8.7 
15.71 



Figure 5: Average values of Shannon entropy, Simpson's index, space-filling ratio and expressive- 
ness with D-stimulaton. 

2.13. Class C 10 

Phenomenology of automata's behaviour is richest amongst all classes. D- 
stimulation leads to formation of trains of excitation wave-fragments, propagating 
in north-west, north-east, south-west and south-east directions; domains of lo- 
calised excitations inside boundaries of D, travelling localizations with and with- 
out traces of refractory states (Fig.|4j5). Typically excited cells are in slight major- 
ity in configurations generated by automata form Cio (Fig. |2} Cio). The class Cio 
includes 61 functions ^{2alb), = 3, . . . , 8; ^{lalb), = 4, . . . , 8. 

2.14. Class C\\ 

Behaviour of automata from class Cn is somewhat similar to that of class C9 
but sometimes localizations travelling along north-south and west-east axis are 
'linked' by segments from excitation wave-fronts (Fig. [4]:). The class C\\ has 
27 functions ^(2312); ^(232a), a = 2,...,8; ^(2412); ^(242a), a = 2,3; 
^{2abc), a = 5, . . . ,S, b = 1,2, b < c <3. 

3. Hierarchies 

Finding 1. Classes of excitable cellular automata with retained refractoriness 
obey the following power hierarchy: 



C4 > Cio > {Ci,Cii} > C7 > C9 > {C3,C8} > Ce > C5 > C2. 



This is a direct consequence of our joining functions into classes. 
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c, c, c, c, c, c, c, c, c, c,„ c„ 

Classes 

Figure 6: Space-filling ratio values for classes of retained refractoriness rules. Classes are indi- 
cated on horizontal axis. Space-filling is shown on vertical axis. All 225 functions are mapped 
to show distribution of space-filling ratio inside each class. There is less discs on the chart than 
functions because some functions from the same class have the same space-filling ratio values. 
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Finding 2. Classes of excitable cellular automata with retained refractoriness 
obey the following hierarchy of space-filling ratio: 

{Ci , C2, C3, C5 , Cg} > Cio > C7 > Ci 1 > {Cg, C9} 

Space-filling ratio values averaged amongst each class are shown in Fig. |5] and 
the distribution of functions inside each class in Fig. [6j Classes C\ to Cg, and Cg 
show small variance in space-filling ratio while class C\\ shows largest variance. 
Roughly the classes can be split into two space-filling groups: high space-filling 
ratio (Ci , . . . , Cg) and low space-filling ratio (C7, . . . , Ci 1). 

Finding 3. Classes of excitable cellular automata with retained refractoriness 
obey the following hierarchy of diversity: 

Ci > C3 > C4 > Cio > C2 > Cg > Cii > C7 > C9 > Cg > C5. 

We can draw an analogy between neighbourhood states in a given config- 
uration of cellular automaton and living species in a population. Thus we can 
measure a diversity of a function by calculating Simpson's index and Shannon en- 
tropy of a population of neighbourhood states in a single configuration generated 
by the function (Fig. |7]). Classes C5 and Cg exhibit the lowest diversity because 
initial random excitation pattern D extinguishes quickly and only single propagat- 
ing wave-front travels away from the site of original excitation. The front leaves 
behind an almost homogeneous domain of refractory states. Classes Ci, C3 and 
C4 hold top ranks in the diversity hierarchy because the initial random patterns ex- 
pands into a growing quasi-random configuration of excited and refractory states. 
In terms of Shannon entropy class C3 shows higher degree of diversity than class 
C4 (Fig. [5]), possibly, because space-time dynamics is biased towards formation 
of extended wave-fronts propagating north, south, east and west, which interact 
with fragmentary wave-patterns travelling north-west, north-east, south-west and 
south-east. Class C\q, which demonstrates high variety of travelling localizations, 
and stays just above the middle of the diversity hierarchy. 

Finding 4. Classes of excitable cellular automata with retained refractoriness 
obey the following hierarchy of expressivenes: 

Cg > C9 > Cii > C7 > Cio > Ci > C3 > C4 > C2 > Cg > C5. 
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Figure 7: Distribution of classes on morphological diversity, (a) Simpson's index vs Shannon 
entropy and (b) Shannon entropy with D-stimulation. 




Figure 8: Space-filling ratio vs Shannon entropy with D-stimulation. 
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Distribution of functions in a space filing versus Shannon entropy space is 
shown in Fig. [8| Despite having below average diversities classes Cg and Cg oc- 
cupy top level of the expressivness hierarchy. This is because their space filling 
demands are incredibly modest. Apart of a random refractory domain left after D- 
perturbation, we observe only propagating excitation wave-front in Cg and trace 
leaving localizations in Cg. They are followed by Cn. Automata from Cn exhibit 
varieties of localised and spatially extended excitations. 

Finding 5. Classes of excitable cellular automata with retained refractoriness 
obey the following hierarchy of generative diversity: 

Ci >C2>C9>Cii >Cio>C7>C4>C8>{C3,C5,C6}. 

Average values of Shannon entropy of configurations generated by point- wise 



excitation are shown in Fig. 1 1 Classes Ci and C2 are at the top of generative 
diversity hierarchy. Typical configurations developed after a point- wise excitation 
include almost-rectangular patterns (Fig. |9^) and rectangular patterns with 2/3 
diamond comers (Fig. |9j5). These patterns are symmetric yet sophisticated in 
terms of visual complexity patterns of excited and refractory states. 

In the descending order of generative diversity, classes C\ and C2 are followed 
by classes Cg and Cn. Patterns generated in the result of point- wise stimulation 
include 

• Sierpinski carpets expanding north and south from longest sides of two-cell 
point-wise excitation (Fig.flOk); 



rhomboid excitation wave fronts, followed by travelling localizations and 
leaving behind a domain of disordered sparsely distributed refractory states 
(Fig.flOb); 



• rhomboid wave-fronts forming ordered strips of refractory states (Fig.[TO|:d); 

• and combinations of localised and spatially-extended refractory-state do- 
mains (Fig.[TO^f). 

Class Cio occupies a mid-range of generative diversity. A point-wise exci- 
tation generates rhomboid patterns of wave packets, each excitation front is fol- 
lowed by a refractory tail (Fig. |9]:d). Some functions of the class also produce a 
narrow band of resting states, surrounded by excitation domain (Fig.[9|:). Classes 
C3, C5, C4, Q Cg have lowest generative diversity. A point- wise excitation leads 
to formation of a circular excitation wave, which leaves a trail of refractory states. 
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(a) ^(1228) e Ci 



(b) ^(1222) e C2 




(c) ^(2413) e Cio (d) ^(2426) G do 



Figure 9: Typical configurations, developed after point-wise excitation, generated by functions 
(a) Ci, (b) C2, (cd) Cio- Names of functions are indicated in the sub-captions. Snapshots of 
configurations are taken 100 time steps after point-wise perturbation. Excited cells are red (light 
grey), refractory cells are blue (dark grey) and resting cells are white. 
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(a) ^(2211) e C9 (b) ^(2218) G C9 (c) ^(2222) e C9 (d) ^(2223) £ C9 




(e) ^(2312) e Cii (f) ^(2322) e Cu 



Figure 10: Typical configurations, developed after point-wise excitation, generated by functions 
(a-d) Cg, and (ef) Cu. Exact functions are indicated in the sub-captions. Snapshots of configu- 
rations are taken 100 time steps after point-wise perturbation. Excited cells are red (light grey), 
refractory cells are blue (dark grey) and resting cells are white. 
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Class 


Shannon entropy 


Simpson's index 


Space-filling ratio 


Expressiveness 


Ci 


6.89 


1.0 


0.90 


7.65 


C2 


5.52 


0.98 


0.84 


6.6 


C3 


0.31 


0.10 


1.0 


0.31 


C4 


2.86 


0.94 


0.08 


34.32 


C5 


0.31 


0.1 


1.0 


0.31 


C6 


0.31 


0.1 


1.0 


0.31 


C7 


2.93 


0.94 


0.09 


32.58 


Cg 


2.71 


0.93 


0.07 


38.79 


C9 


4.55 


0.96 


0.22 


20.59 




3.16 


0.94 


0.27 


11.71 


Cii 


3.80 


0.95 


0.21 


17.91 



Figure 11: Average values of Shannon entropy, Simpson's index, space-filling ratio and expres- 
siveness calculated on configurations generated by a point-wise excitation. 



Space filling vs. Shannon entropy (point-wise excitation) 
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Figure 12: Space-filling ratio vs Shannon entropy. Point-wise excitation. 
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(a) (b) (c) (d) (e) (f) (g) (h) 



Figure 13: Examples of travelling localizations found in configurations of automata with retained 
refractoriness, (a-e) Gliders, (f-h) Trains. Excited cells are red (light grey), refractory cells are 
blue (dark grey) and resting cells are white. 



Finding 6. Classes of excitable cellular automata with retained refractoriness 
obey the following hierarchy of generative expressiveness: 

C8>C4>C7>C9>Cii >Cio>Ci >C2>{C3,C5,C6}. 



See detailed structure of the hierarchy in Fig. 12 



4. Travelling localisations 

Travelling localised excitations is a unique phenomenon of sub-excitable me- 
dia, widely used in theoretical and experimental laboratory implementations of 
novel computing architectures, see e.g. Most common localizations found 
in experiments with excitable automata with retained refractoriness are gliders, 
trains, glider guns, and spaceships. These structures are well-known in Conway's 
Game-of-Life cellular automata (and we are using terminology of the Conway's 
Game of Life) however they display a range of uncommon properties when inhab- 
iting discrete excitable medium. 

Examples of gliders (localizations which do not leave traces) and trains (lo- 
calizations which leave traces made of stationary refractory patterns) are shown in 
Fig. |4| all localizations shown there travel north. A glider shown in Fig. con- 
sists of three excited and states. The small gliders can be joined together to make 
propagating wave-fragments: in Fig.|4|5 we see a glider made of three excited and 
three refractory states joined with a localisation of two excited and two refractory 
states. The gliders move along columns and rows of a cellular array. Similarly 
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to our previous models of sub-excitable discrete media, e.g. we found only 
one excited localisation which travels along diagonals of the cellular array. The 
diagonal glide consists of three excited and three refractory states, an exemplar 
snapshot is shown in Fig. [4]:. 

Gliders which leave traces of stationary patterns behind are called trains. The 
trains are very common in classes Cg, Cio and Cn. A minimal train consists of 
a pair of excited states leaving a column/row of refractory pairs behind (Fig.|4|l). 
In terms of trace minimality, the train shown in Fig.|4^ is minimal localisation of 
three excited states which leaves a trace of refractory states one-cell wide. Trains 
can be linked with gliders and other trains, an example is shown in (Fig.|4]r). Trains 
leaving localised traces of refractory states are shown in (Fig. |4^h). 

Travelling localizations which generate streams of other localizations, e.g. 
gliders, are called mobile guns (Fig. 14). A mobile gun can emit localisation 
in any directions, subject to the lattice's discreteness, apart of the direction of the 
gun's motion. Here we illustrate just few examples: gun travels east and emits 
trains travelling south (Fig. [14k), gun travels east and emits gliders travelling west 
(Fig. 14 d), gun travels north and emits trains travelling east (Fig. 14 gun travels 
north and emits gliders travelling west (Fig. [T4|i). A complex gun is shown in 
(Fig. 14 j): the gun travels north and emits three streams of gliders, one stream 
travels west and two streams travel east. 

Interactions between gliders is pretty rich, gliders and trains can reflect or 
annihilate in the result of the collision, also refractory stationary patterns can be 



formed in the results of such collisions. An example is shown in Fig. 15 Two 



mobile guns travel north. They both leave a trace of refractory states, as trains do. 
Westward gun emits a single stream of gliders travelling east. Eastward gun emits 
two streams of gliders. Eastward stream of westward gun collides with westward 
stream of eastward gun. Refractory patters are formed in the result of collision, 
and glider streams are pruned: only every fifth glider of westward stream survives 
collision and only every fourth glider of eastward stream survives the collision 
(Fig.[T5]). 

Functions ^{2alb), where a = 3,..., 8 and b = 1,2, generate spaceships 
(b=l) and Sierpinski flying carpets (b = 2). Sierpinski carpets and spaceships ob- 
served in many types of cellular automata including Conway's Game of Life lfT9l 
[16]. However, we believe, travelling Sierpinski carpets have never been observed 
before in cellular automata. Examples of Sierpinski carpets generated by function 
^(231 1) from seeds of different sizes are shown in Fig. 16 The carpets always 
leave a trail of refractory states. Therefore they can be considered as analogies of 
puffer trains fT5] in the Game of Life cellular automata. 
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Figure 14: Examples of mobile glider guns found in development of automata with retained re- 
fractoriness. Excited cells are red (light grey), refractory cells are blue (dark grey) and resting 
cells are white. 
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Figure 15: Example of interactions between travelling excitations in automata with retained re- 
fractoriness. Excited cells are red (light grey), refractory cells are blue (dark grey) and resting 
cells are white. 
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Figure 16: 'Flying' Sierpinski carpets, generated by automata of class C9. Carpets travel north. 
Patterns (from the left to the right) are generated by lines of two, three, four and five excited cells. 
Launching pad is constructed by a glider, which travels east and leaves a refractory trail. The 
configurations are generated by function ^(2311). Excited cells are red (light grey), refractory 
cells are blue (dark grey). 
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(a) 





(c) 







id) 

Figure 17: Formation and propagation of spaceships in automata of class Cn. Patterns (from the 
left to the right) are generated by lines of two, three, four and five excited cells. Launching pad is 
constructed by glider travelling east and leaving refractory trail, (a) Generation stage, (be) Take- 
off stage, (d) Flying. The configurations are generated by function ^(2312). Excited cells are 
red (light grey), refractory cells are blue (dark grey). 
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Spaceships do not leave refractory debris behind. Several stages of space- 



ships development are illustrated in Fig. 17 A strip of excited states, e.g. '+ 



"++', "+++' etc., leads to formation of self-similar growing rhomb (Fig. 17 1 and 
Fig. 10 le). When the south part of the rhomb gets in touch with a refractory 
strip it becomes 'infected' with refractory states, or frozen, thus only north part of 
the rhomb continues travelling (Fig. [ITj^c). At some stage a spaceship becomes 
disconnected from refractory debris it left and gets 'airborne' (Fig.[T7|j). 

Both Sierpinski flying carpets and spaceships observed are rather close rela- 
tives of 'grayships' [16J, because they are comprised of still patterns which can 
be extended to arbitrary size. The size of each carpet or spaceship is determined 
by the distance of its seed from the launchpad (a chain of refractory states, which 
prevents the carpet/spaceship from expanding in one direction). 

Collisions between the flying carpets and spaceships develop along 'classical' 
scenarios of collisions between travelling localizations: annihilation and reflec- 
tion, possibly with emergence of new localizations (Fig.fTSl). 



5. Discussion 

Aiming to design a cellular automaton model of a discrete memristive medium 
we developed a model of excitable cellular automaton with retained refractoriness. 
The model is ideologically similar to our previous model of cellular automata with 
retained excitation [i5j| but the focus shifted to retained refractoriness. An excited 
cell takes refractory state unconditionally, i.e. independently of the number of its 
excited neighbours, however the refractory cell can remain in its refractory state 
depending on the number of its excited neighbours. 

A propagation of excitation in a medium can be considered as an analog of 
electrical current. A refractory state is an automaton analog of high-resistance 
state of a memristive element, and a switching form the refractory state to rest- 
ing state, controlled by a density of local excitation, is an analog of switching 
the memristive element to its state of low-resistance. The automata studied are 
semi-memristive because switching from low-resistance (resting) state to high- 
resistance (refractory) state is not controlled by electrical current (density of local 
excitation). 

Classical automaton models of excitable media are usually based on the prin- 
ciple of threshold excitation: a resting cell excites if the number of excited neigh- 
bours exceeds certain threshold. This threshold-view on excitation dynamics was 
the main reason why mobile localizations have not been discovered for so long 
time. In [1] we introduced interval-based excitation: resting cell takes excited 
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(a) (b) (c) (d) 




Figure 18: Interaction of flying Sierpinski carpets in automata of class Cn. (a) Carpets approach 
each other, (be) Even distance: (b) full body impact and (c) formation of debris/annihilation, 
(de) Odd distance: (d) full body impact and (e) formation of new travelling localizations. The 
configurations are generated by function ^(2811); seed's size is two cells; distance of a seed 
from southern refractory strip (launchpad) is 50 cells. Excited cells are red (light grey), refractory 
cells are blue (dark grey). 
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state if the number of its excited neighbours belongs to some specified interval. 
The approach proved to be productive — a range of non-classical computing de- 
vices [4J, based on the rich phenomenology of excitation in interval-excitation- 
based models, was designed. Therefore we employed interval-based excitation 
and also interval-based recovery, i.e. transition from refractory to resting state, to 
imitate semi-memristive automata. 

Based on functions' behavioural grouping, visual categorisation of configura- 
tions generated in a response to random perturbation, ordering functions on their 
morphological diversity and expressiveness we selected eleven classes of func- 
tions: Ci, . . . ,Cii. 

The functions of classes show high degree of space-filling ratios, 

while classes C7, . . . , Ci i mid-range and low space-filling ratios. Low space-filling 
classes can be further sub-divided into classes with extinct excitation, namely C7 
and Cg, where only a single centripetal wave-front is formed, and classes that sup- 
port travelling localizations: Cg.Cio.Cn. The low space-filling classes C7, . . . ,Cii 
also demonstrate higher ratio of excited states over refractory states in configura- 
tion functions of the classes generate. 

Top four classes in the hierarchy of morphological diversity are Ci, C3, C4 and 
Cio. Automata from Ci, C3 and C4 show rather quasi-chaotic response to spatially 
extended random excitation, and this is the reason for their high morphological 
diversity. Class Cio is more interesting. A spatially extended random perturba- 
tion leaves a random configuration of refractory states filled with non-expanding, 
sometimes breathing, domains of localised excitations and a combination of prop- 
agating wave-fronts, wave-fragments and travelling localizations. 

Top four classes in the hierarchy of generative diversity are Ci, C2, C9 and Ci 1. 
Classes Ci and C2 show quasi-chaotic dynamics and even a point-wise excitation 
of a resting lattice leads to formation of a complex growing pattern. Class Cg got 
to the top four of generative diversity because a point-wise excitations give birth 
to expanding Sierpinski carpets, and disordered and highly-ordered arrangements 
of refractory states. 

To refine our classification even further we introduced a measure of expres- 
siveness, which relates morphological diversity to space-filling ratio, and thus re- 
flects an 'economy of diversity' in cellular automaton configurations. As you can 
see in Fig. [8} classes C9, . . . Cn show higher degrees of expressiveness. Thus we 
can propose that functions which support travelling localizations may be at the 
mid-range levels of morphological diversity but top levels of expressiveness. 

We briefly characterised a range of travelling localizations which emerged 
in the cellular automata studied. We demonstrated that many localizations are 
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capable for laying sophisticated patterns of refractory, high-resistive in memris- 
tor's terminology, states. Such patterns could be used in further studies to com- 
partmentalise otherwise homogeneous spatially extended memristive media. We 
found outcomes of collisions between propagating localisation are quite rich and 
therefore all scenarios of collision-based computing [3| can be implemented in 
semi-memristive cellular automata. Our further studies will be concerned with 
developing fully memristive cellular automaton model and designing laboratory 
experimental prototypes of large-scale memristive media. 
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